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NOMENCLATURE 

E 

modulus  of  elasticity  of  beam,  lbf/in. 

I 

moment  of  inertia  of  beam  for  axis  perpendicular  to  the 
forcing  plane,  in. 

Ab 

2 
cross  sectional  area  of  beam,  in. 

L 

length  of  beam,  in. 

P 

density  of  beam  material,  lbm/in. 

Fi 

amplitude  of  the  7.  component  of  the  beam  forcing  function, 

lbf/in. 

F 

beam  forcing  function,  lbf/in. 

7 

circular  frequency  of  vibration,  rad/sec. 

f 

frequency  of  vibration,  cycles/sec. 

m. 

generic  pendulum  mass,  lbm. 

Jl 

generic  pendulum  length,  in. 

X 

beam  coordinate  from  fixed  end,  in. 

v.(x,t) 

beam  displacement,  in. 

v  '(t) 

1 

displacement  of  i   generic  pendulum  mass,  in. 

t 

time,  sec. 

6 

acceleration  of  gravity  in. /sec. 

V 

shear  in  beam   lbf . 

M 

bending  moment  in  beam    in.  lbf. 
n 
)  m.   ,   lbm 

in 

P 

J=1 

^V2           -i                                                                 .'       1 

V  EI    '   in* 

iv 

PART  I 


INTRODUCTION 

When  tall,  flexible  structures  such  as  light  poles,  antenna  masts;  and 
stacks  are  subjected  to  transverse  vibratory  loadings,  large,  undesirable 
deflections  sometimes  result.  Such  vibratory  loadings  are  most  often  caused  by 
winds  of  high  speed  which  cause  a  vortex  shedding  effect.  A  vortex  shedding 
effect  of  this  nature  causes  a  dynamic  loading  pattern  on  the  structure  in  the 
cross-wind  direction.  At  the  present  time,  little  is  known  about  the  nature 
of  these  forces  except  that  it  depends  somewhat  on  wind  speed  and  structural 
geometry. 

Reed  and  Duncan  (ref .  l)  have  studied  methods  to  suppress  wind  induced 
oscillations  of  the  nature  just  described.   Such  structures  as  the  Saturn  V 
launch  vehicle  and  Navy  antenna  systems  were  considered.   In  both  instances, 
oscillations  were  suppressed  by  attaching  energy  absorbinb  elements  to  the 
structure. 

In  the  case  of  the  Saturn  V  launch  vehicle,  dampers  were  placed  between 
the  launch  vehicle  and  the  tower.  Results  of  a  wind  tunnel  study  of  a  scaled 
model  revealed  that  the  dynamic  response  of  the  vehicle  was  relatively  low 
even  at  wind  speeds  of  60  knots . 

In  the  case  of  the  Navy  antenna  systems,  a  cluster  of  plastic  coated 


"Wilmer  H.  Reed  III  and  Rodney  L.  Duncan,  Proceedings  of  the  Tenth 
Midwestern  Mechanics  Conference.   (Fort  Collins,  Colorado,  I967),  pp.  88l 
897. 


chains  surrounded  by  a  neoprene  shroud  was  attached  at  the  upper  end  of 
the  antenna.  The  cluster  was  placed  inside  the  cylindrical  antenna  and 
served  as  an  effective  impact  type  damper.   In  both  of  these  examples, 
the  vibration  absorbers  were  designed  to  absorb  energy  so  that  large 
oscillations  would  not  occur. 

The  purpose  of  this  paper  was  to  study  the  effects  of  a  pendulum 
type  vibration  absorber  which,  when  properly  designed,  would  make  the 
absorbed  system  insensitive  to  certain  forcing  frequencies.  As  in  the 
case  of  the  previous  examples,  it  was  felt  that  this  type  of  vibration 
absorber  would  be  well  suited  for  tall  flexible  structures  subjected  to 
dynamic  loads  in  the  transverse  direction.  The  purpose  of  the  absorber 
was  not  one  of  energy  absorbtion.   Instead,  the  pendulum  absorber  was  to 
change  the  resonant  frequencies  of  the  absorbed  system  so  that  they  would 
not  correspond  to  the  forcing  frequency  or  frequencies.  Both  single  and 
multiple  pendulum  absorbers  were  considered. 

Since  internal  damping  in  structures  is  usually  small,  it  was 
neglected  in  this  paper.   It  was  also  assumed  that  the  pendulum  absorb- 
ers were  frictionless,  and  that  displacements  were  sufficiently  small  so 
that  a  linear  analysis  could  be  used. 

When  multiple  pendulums  with  flexible  links  are  excited  above  the 
fundamental  mode,  the  displacements  tend  to  become  nonplaner.  It  was 
assumed  in  this  paper  that  the  pendulum  displacements  were  restricted  to 
the  plane  of  the  forcing  function.  This  could  be  done  by  suitable  link 
design.   It  was  also  assumed  that  the  masses  of  the  pendulums  could  be 
considered  as  point  masses.  The  structure  which  was  analyzed  in  this 
paper  was  a  vertical,  Hookean,  isotropic  homogeneous  cantilever  beam  of 


constant  cross  section.  The  absorber  was  placed  at  the  unfixed  end.  It 
was  also  assumed  that  the  forcing  function  was  uniform  along  the  length 
of  the  beam,  and  that  it  was  composed  of  a  finite  number  of  sinusoidal 
components.  This  model  is  illustrated  in  figure  1  of  chapter  2. 

In  this  paper,  displacement  equations  were  developed  for  various 
absorbed  cantilever  systems.  From  these  displacement  equations,  charac- 
teristic equations  (equations  whose  roots  yield  resonant  frequencies) 
were  developed.  Computer  programs  written  in  Fortran  h,   Level  G,  were 
developed  to  evaluate  these  displacement  equations  and  determine  the 
resonant  frequencies  of  the  various  cantilever  systems. 

It  was  discovered  that  if  none  of  the  forcing  function  frequencies 
are  resonant  frequencies  of  the  fixed-pinned  beam,  then  it  is  possible 
to  design  an  absorber  such  that  the  absorbed  cantilever  beam  will  be 
effectively  pinned  at  the  absorbed  end.   It  can  also  be  guaranteed  that 
resonance  will  not  occur  at  these  frequencies.  These  statements  are 
proved  in  Chapter  3- 

Experiments  were  performed  in  order  to  verify  the  validity  of  the 
characteristic  equations.  Chapter  h   contains  the  description  and  results 
of  these  experiments. 


PART  II 
DEVELOPMENT  OF  DISPLACEMENT  EQUATIONS 

In  this  chapter,  the  displacement  equations  for  various  cantilever 
beam  systems  are  developed.   In  particular,  displacement  equations  for  a 
cantilever  beam  with  a  free  end,  with  a  simple  pendulum  absorber,  and 
with  a  multiple  pendulum  absorber  are  derived.   In  each  case,  it  is 
assumed  that  the  beam  experiences  a  forcing  function  which  is  uniform 
along  its  length  and  which  consists  of  a  finite  number  of  sinusoidal 
components.  This  situation  is  illustrated  in  figure  1. 

A  free  body  diagram  of  a  differential  element  of  the  beam  at  a 

point  x  is  shown  below  in  figure  2.  Equating  the  sum  of  the  forces  in 

the  "v"  direction  to  the  product  of  the  mass  of  the  differential  element 

and  the  acceleration  of  the  differential  element  in  the  "v"  direction, 

the  following  equation  is  obtained: 

N2 
j 

oV 


Fdx  +  V  -  (v  +  g-  dx)  -  PAb  &if&   dx  .  (2-1) 


V  -  EI  ^feJil  #  (2.2) 

dx5 

Using  equation  (2-2),  the  equation  from  elementary  beam  theory  which 

relates  shear  and  displacement,  equation  (2-1 )  reduces  to 

auv(x,t)  .  ^k  a2v(x.t)  _  jl  ,p     , 

dx4    +  EI     at2    ~  EI   *  ^2"5a) 


The  boundary  conditions  for  this  problem  are  as  follows : 
1.  v(0,t)  -  0  . 


(2->b) 


2. 


3- 


dv(x.t) 


dx 


=  0  . 


x=0 


o  v(x.t) 
^2 


=  0  . 


x=L 


k.     EI 


aVx.t) 


^ 


=  V 


x=L 


x=L 


(2-3c) 


(2-3d) 


(2-3e) 


In  appendix  A,  the  principle  of  superposition  is  proven  for  this 
problem.  This  principle  states  that  if  v.(x,t)  is  the  displacement  function 
for  a  beam  acted  upon  by  a  forcing  function,  F.  sin7.t,  then  the  resultant 
displacement  function  of  a  beam  acted  upon  by  a  forcing  function 


Z    F. sin7st,  is  as  follows: 
i=l      1 


r 
r(x,t)  =   )   Vi(x,t)   . 


(2-M 


i=l 


Due  to  this  principle,  the  derivations  of  the  displacement  equations  will  be 
made  by  assuming  a  forcing  function  of  a  single  sinusoidal  component 

F(t)  =  Fsin7t  . 

Equation  (2 -3a)  now  takes  the  following  form: 


A(x.t)  ,  °\    b2v(x.t)       _F_ 


EI 


EI 


sinyt  . 


dx      **    dtc 
By  assuming  v(x,t)  =  r(x)p(t),  equation  (2-5)  reduces  to: 


(2-5) 


P(t)  &£i  +  ^M  ism. .  i  8inn  .  (2.6) 

dx  dt 

The  solution  of  equation  (2-6)  is  the  sum  of  a  homogeneous  solution  and  a 
particular  solution.  Due  to  a  slight  amount  of  damping,  which  is  present 
in  all  physical  systems,  the  steady  state  solution  will  contain  only  the 
particular  solution. 

For  the  steady  state  solution,  it  is  assumed  that  p(t)  is  defined  by 
equation  (2-7). 

p(t)  =  CjSinrt  +  C2cos7t  .  (2-7) 

By  combining  equations  (2-6)  and  (2-7),  the  following  equation  is  obtained: 

,       2 
(Clsinrt  +  C2cos7t)  [LE&l  -  Z^k  r(x)]  =  WjBk   §  (g^} 

Since  the  right  side  of  equation  (2-8)  is  a  function  of  t  only;  the  term  in 
brackets  must  be  a  constant,  and  Cp  must  be  zero.  By  letting 


-  dx 


-  Q  , 


C-^  must  take  the  form  expressed  in  equation  (2-9)  . 

Cl  =  «H  '  (2-9) 

The  solution  for  r(x)  now  becomes: 

r(x)  =  Asinpx  +  Bcospx  +  Csint^x  +  Dcoshpx  -  ~%  .  (2-10) 

2  P 

P     =     EI       ' 

Since  A,  B,  C,  D,  and  Q  are  constants  which  are  yet  to  be  determined,  no  loss 

of  generality  will  result  by  expressing  v(x,t)  as  in  equation  (2-ll). 


v(x,t)  =  sin7t  J  AsinPx  +  Bcospx  +  CsinhPx  +  Dcoshpx  ~r \  .     (2-11 ) 

L  fTEI  J 

By  substituting  this  expression  in  equations  (2-3b),  (2-3c),  (2-3d), 

and  (2~3e)j  the  four  boundary  condition  equations  in  terms  of  A,  B,  C,  and 

Dare  shown  in  equations  (2-12a),  (2-12b),  (2-12c),  and  (2-12d). 

B  +  D  -  -]~  =  0  .  (2-12a) 

P  EI 

A  +  C  =  0  .  (2-12b) 

.    -AsinpL  -  BcospL  +  CsinhPL  +  DcoshPL  =  0  .  (2-12c) 

EI05  -NacosPL  +  BsinpL  +  CcoshpL  +  DsinhpLf     =  V  .  (2-12d) 

I  ■>  |x=L 

A  simultaneous  solution  of  these  four  equations  yields  equations   (2-13a), 
(2-13b),    (2-13c),  and  (2-13d). 

_(v|        )   (cosfSL  f  coshpL)  +  _E_|coshpLsinpL  +  C0S3LsinhPL} 

A  ..  x  L        El34sinyt  p\l  I i  /2.„   } 

A  ~  2  +  2coshpLcos3L  '         Kd  ±5a; 


-£—  cospL  -  A(sinPL  +  sinhpL) 

B  =  01 cosPL  +  coshPL '  <2"15b) 

C  =  -A   .  (2-13c) 

D  =  -f—    -  B  .  (2-13d) 

0  EI 

A  convenient  form  for  the  displacement  equation  can  be  made  by  substi- 
tuting equations  (2-13a),  (2-13b),  (2-13c),  and  (2-13d)  into  equation  (2-ll), 
and  by  expressing  the  displacement  equation  in  terms  of  A.  In  this  manner, 
various  displacement  equations  which  arise  from  the  consideration  of  different 


absorbers  may  easily  be  formed. 

By  performing  these  steps,  the  displacement  equation  becomes  equa- 
tion (2-13e). 


v(x,t)  =  sin7t 


A(B  )  +  -f~  (Bp) 

P  EI     • 


(2-13e) 


B^sinPx  -  sinhPx  -  (gg  + "  gggg  (cosPx  -  coshPx)      (2-l5f  ) 

^    cosPL  coshPx  +  cosPL  coshPx   ,  ,~  ..,  \ 

B2  =      cosPL  +  coshPL        "  1   '  (2"15g) 

If  the  beam  is  free  at  the  end,  V|  _T  =0;  and,  by  using 

equations  (2-13),  the  displacement  equation  becomes  equation  (2-l4) . 


Fsinyt 


v(x,t)  - 

P  EI 


(GjHGj,)  +  C5)J  .  (2 -lta] 

_    JcoshpL  sinPL  +  cosPL  sinhPLL  (o   -,iK\ 

1~\  2  +  2  cosPL  coshpL    J  *  ^-l»l>; 

CQ  =  sinpx  -  sinhPx-  (§^%~r3JJ^)( e°spx  -  coshPxJ  (2-l^c) 

2  \cosPL  +  coshPL/  \  / 

c  .  JcoshpL  cosPx  +  cosPL  coshpx  .  |  (   ^  j 

3  ^     cosPL  +  coshPL  J  v     ' 

By  studying  equations  (2-lUc)  and  2-lUd),  it  becomes  apparent  that  both 
terms  Cp  and  C,  are  bounded  when  P  and  x  are  positive  and  bounded.  It  is 

therefore  evident  that  for  resonance,  C,  must  become  infinite.  Resonance 

for  a  free-ended  cantilever  will  therefore  occur  when  equation  (2-15)  is 

satisfied. 

2  +  2  cosPL  coshpL  =   0  .  (2-15  ) 

Equation  15  is  transcendental  and  hence,  values  of  P  cannot  be  solved  for 

directly.  A  computer  program  was  used  to  find  the  first  15  roots  of  Eq. 
(2-15). 


10 


x=L 


These  results  appear  in  appendix  C. 

If  a  simple  pendulum  is  attached  to  the  upper  end  of  the  column,  V| 
is  no  longer  zero,  and  a  different  displacement  equation  will  result. 
Consider  a  pendulum  of  length  I     and  mass  iru  .   Equation  (2-l6)  is  obtained 
from  a  free  body  diagram  of  m  .   (See  figure  }) 


v(L,t) 


m±B    ¥ 


Figure  3« 
Free -Body  Diagram  of  Single  Mass  Pendulum 


m 


?  yi^     mig 


1    at2 


v(L,t)  -  v±   (t) 


(2-16) 


Equation  (2-l6)  simplifies  to  equation  (2-17)  when  r(L)sin7t  is  substi- 
tuted for  v(x,t). 


Bvi(t)    K     ' 

at2       li    x 


f-  r(L)sin7t 


(2-17) 


The  solution  of  this  equation  consists  of  the  sum  of  a  homogeneous  solution  and 


11 


a  particular  solution.  Due  to  a  small  amount  of  damping,  only  the  particular 

solution  will  remain  in  the  steady  state  solution.  By  assuming  the  form  of 

v'(t)  expressed  by  equation  (2-l8),  the  steady  state  solution  for  v'(t)  is 

found  to  be  as  expressed  in  equation  (2-19). 


v|(t)  =  C^sinrt  +  C2cos7t 


(2-18) 


1  tj.  \      r(L)sinyt 


(2-19) 


The  shear  force  V|  _T  at  the  upper  end  of  the  beam  must  be  equal  to 


'x=L 


in. 


3t2 


Thus, 


-m,r(L)sin  t 


1  - 


V 


2    =VU  • 


(2-20) 


G 


After  substituting  equation  (2-20)  into  equations  (2-13 )  and  simplifying  the 
results,  the  displacement  equations  for  the  beam,  and  the  pendulum  mass  reduce 
jto  equations  (2-2l). 


v(x,t)  =  £-p2&  (Di)(l)2)  +  ^ 
P  EI 


ij  ' 


(2-21a) 


rUt)   -i%^  . 


(2 -21b) 


12 

A   2  coshpL  cospL  "\  (     "ml   N 

(2 -21c) 

^11   V    cos^L  +  coshpL  )\          l1J     ' 

72    g 

EI35(coshPL  sinpL  +  cos&L  sinhPL) 
12  "          cosPL  •»•  coshpL 

(2-21d) 

^ 

(2-21e) 

«  "  i   ■  i  ' 

2    g 

^ .  [*•  -  **«  -  (£g :  &g)  (-*—*)  ] 

(2-21f) 

• 

EI35(2  +  2  "cosh0L  cos(3L) 
15  '      cospL  +  coshpL 

(2~21g) 

D       Dll  +  D12 

(2-21h) 

1  '  V  V  +  Dl5  * 

D2  -  ["si^x  -  sinhfa  -  (gg;:;^)(«o^  -  ooshPx)  ]  .        (2-2U) 


coshftL  cosPx  +  cosftL  coshpx 
cosPL  +  coshPL 


(2-21J) 


Since  cospL  +  coshPL  £   0  for  BL  >  0,  D,  will  remain  finite  for  finite 
values  of  x.  Consequently,  if  resonance  occurs,  the  product  (D,  )(Dp)  must 
grow  without  bound.  Since  D^  in  equation  (2-2l)  is  bounded  for  positive 
finite  values  of  0  and  x,  resonance  must  occur  when  D.  grows  without  bound. 
By  observing  the  equations  which  define  D..,  it  is  evident  that  the  characteris- 
tic equation  for  a  system  consisting  of  a  cantilever  beam  and  a  simple  pendu- 
lum attached  at  the  end,  must  be  equation  (2-22). 


15 


0  =  D 


1   i 
,2    6 


Ik 


m. 


D 


15  * 


(2-22) 


Equation  (2_22)  is  transcendental  and  hence  the  roots  cannot  be  obtained 
directly.   It  is  interesting  to  note  that  not  only  the  beam  parameters,  but 
also  the  pendulum  parameters  influence  the  roots  of  equation  (2-22),  and  hence 
influence  the  spectrum  of  resonant  frequencies  of  the  coupled  system.  A 
computer  program  was  used  to  solve  for  the  roots  of  equation  (2-22).  A  list- 
ing of  the  program  is  contained  in  appendix  D. 

If  a  multiple  pendulum  as  shown  in  figure  k   is  attached  to  the  free  end 
of  a  cantilever  beam,  the  shear,  V,  at  x=L  must  be  as  expressed  by  equation 
(2-2J).  The  value  of  V  at  x=L 


Cm 


6mr 


\ 
\ 
\ 
V 

\ 


0  m 

n 
0  m 


n-1 


n 


Y(t) 


Figure  h. 
Diagram  of  Beam  and  Multiple  Pendulum  Absorber 


Ik 


'X-L, 


a2Yj(t) 


1=1 


mi     2 
1  btd 


(2-23) 


will  depend  upon  the  vibration  of  the  multiple  pendulum.  Since  the  beam 
displacement  equation  is  dependent  upon  V|   ,  the  attachment  of  such  a 
pendulum  will  completely  alter  the  vibration  characteristics  of  the  beam. 

In  analyzing  the  displacement  equation,  it  is  convenient  to  define 
the  following  notation: 

v. '(t)  =  displacement  of  m.   .  (2-2I4-) 


m. 


in 


n 

=  Z  \  • 


(2-25) 


k=i 


The  parameter  m.   defined  in  equation  (2-25)  is  convenient  for  expressing 
the  tensile  forces  in  the  lengths  of  the  pendulum.  The  tensile  force 
T.  in  link  i  now  becomes  as  expressed  in  equation  (2-26). 


T.  =  g  m. 
1   B  In 


(2-26) 


For  mass  m,,  a  free-body  diagram  yields  the  following: 

v2'(t)  -  v1'(t> 


d2v, ' (t ) 

mi  ~zr~  -  \ 


m 


at£ 

o\'(t) 

"17"   :g1vi(t) 


■v(L,t)  -  vx  (t) 
_ 


+  T„ 


.(2-27) 


ra-i    nu,  ' 
In    2n 

*1     "V 


+  v2  (t) 


"m 


2n 


2  J 


mlnv(L,t) 


(2-28) 

For  a  generic  mass  m.,  which  is  an  internal  mass,  a  free-body  diagram 

yields: 

■v^ft)  -  v'(t) 


m.- 


\  (t)      rvx(t)  •  vi(t) 


1  at 


_  _  m 
i 


+  T.  . 
l+l 


i+1 


]• 


(2-29) 
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d  v!(tj     r  rm.  n 

^~^^H-i(t)LirJ+vi(t) 


min  ,Vn 

V     *i+l  - 


-W^IJ 


feH) 


i+1 
(2-30) 


For  the  end  mass,  m  ,  a  free-body  diagram  yields: 


m 


d2v'(t) 


n 


^' 


G 


fi(t) 


-m 


n 


+  v 


■.-!<*> 


m 
_n 

i 

■-  n  J 


(2-31) 


The  displacement  equations  of  the  "n"  mass  pendulum  become  more  manage- 
able by  using  matrix  notation.  By  using  matrix  notation,  equations  (2-28), 
(2-30),  and  (2-31)  reduce  to  the  following  matrix  equation: 


(M)  (v'(t)}  =  g(N){v.[(t)}  +  {T}  sinn 


»- 


r(L)raln  g 


(2-32) 
(2.32a) 


(M)  - 


m.,  0,  .  .  . ,  0 
0  ,  nu,  0,..,  0 


(2-33) 


0  , .  .  .  ,  0,  m 


n 


16 


(N)  - 


In    2n 


o,  ■••  ,  0, 


m 


2n 


ia  ■  ia 


0 


!^in  ,  ^in 


.   ,  0 


"mi+l  n  ,  mi+l  n   0, '".  0 

£■   i    l 

i+1 


i+1 


(2-3«0 


0,    •••  ,  0, 


m 


-m 


nn 


nn 

n       n 


The  solution  of  equation  (2-32)  consists  of  the  sum  of  the  homoge- 
neous solution  and  the  particular  solution.  Due  to  a  small  amount  of 
damping  in  the  pendulum,  the  steady  state  solution  will  only  contain 
the  particular  solution.  By  assuming  for  the  particular  solution  a 
form  for  (v.  (t))  as  given  in  equation  (2-35),  the  matrix  differential 
equation  is  reduced  to  equation  (2-36). 

{v*(t)}  -  {u}  sin7t  .  (2-35) 


(-/(M)  -  g(N))  (u)  =  {T}  (2-36) 

By  making  use  of  equations  (2-37)  and  (2-38),  equation  (2-36)  be- 
comes equivalent  to  equation  (2-39)- 


(R)  -  (-/!(M)  -  g(N) 

(S)  =  (R)"1  .  ** 
(u)  =  (S)  {T}  . 


(2-37) 

(2-38) 
(2-39) 


By  denoting  the  elements  of  matrix  (S)  by  S.  .    ,    equation  (2-39)  becomes 


(u)   -  I&lp*UL  (s.i} 


(2-1*6) 


*  Matrix(N)  contains  a  diagonal  band  which  is  three  elements  wide.  The 

first  and  last  row  have  only  two  elements  as  indicated. 
**The  case  when  (R)  is  singular  will  be  treated  later. 
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In  equation  (2-40),  (S..)  is  the  column  vector  which  has  elements  identical 
with  the  first  column  of  the  square  matrix  (s). 

By   considering  equations  (2-23),  (2~35),  and  (2-1+0),  it  is  evident  that 
the  expression  for  V|  .  becomes  equation  (2-42). 


V|X=L  "  -?2sin7t 


n 

y  n^  .  (2-^1) 


i=l 


2  ,T  »     sin7t    n 
-7  r(L)m.  g       c- 


y  _lS 


x=L  -       JL         I  Vil  '  (2^2) 
X        i=l 

From  equations  (2-13),  the  displacement  equation  for  a  cantilever  beam  with  a 

multiple  pendulum  (see  figure  h)   becomes  equation  (2-43). 

v(x,t)  -  S-fiart  k  Bl  +  b2]     .  (2-1.38) 

h  -  ••**  -  •**«  -  (as :  ago  Gospx  - coshpx)  •  (2^3b) 

■r         eosftL  coshftx  +  cosPL  coshfix       .,  (p  h*    \ 

*2  =  cos3L  +  coshpL  "  "  X     '  ST***! 


■v2  n 

7  m4„  g 

All  =  " 


^     ym.s.1(-2cgfPLco^+l}     .  (2-43d) 

..  /j   ill  [  cosPL  +  coshpL        J  v 

1    i=l 


l12   cosPL  +  cosh3L 


"+  co  hSL  jcosBL  sin?L  +  cosPL  sinhpLf  .  (2-43e) 


7  m 


A21 


i=l 

(2-43f) 


• 

_  _ 

18 

EI35(2  +  2  coshpL  cosPL) 
22  "     cosPL  +  cosh^L 

(2-ltfg) 

All  +  A12 
1   A21  +  A22 

(2-ltfh) 

By  observing  equations  (2 

-1*3),  it  is  apparent  that  B,  and  Bp 

are 

both 

bounded  for  values  of  x,  P,  and 

L  which  are  positive  and  bounded. 

Therefore, 

if 

resonance  occurs,  then  A,  must  grow  without  bound. 

Since  the  t 

srm 

>  m, s  n  grows  without  bound 

for  certain  frequencies; 

it  is  best 

to 

consider 

I— 

i= 

=1 

as 

the  characteristic  equation 

of  the  system,  equation 

(2 -Ma). 

0  ■ 

"  Dl  +  D2  ' 

(2-kka.) 

h  ■  {slnPL  -  .IrtgL  -  (ffg  t  ggjg)    (cosPL  -  cosW)}       .  (*-*» 


-ijEI^(2  +  2  coshpL  cospL)) 
D2.-J _  (8-1*. ) 


V  mln  /   miSil/  (  C0SPL  *  coshPL  ) 


i=l 

This  chapter  dealt  only  with  the  derivations  of  the  displacement  equa- 
tions and  the  characteristic  equations  for  the  various  cantilever  systems 
under  study.  In  chapter  3>  these  equations  will  be  used  to  study  the  effects 
of  pendulum  parameters  on  the  vibration  characteristics  of  the  cantilever 
systems . 


PART  III 

- 

METHOD  OF  ABSORBER  DESIGN 

Suppose  that  a  certain  fixed-free  beam  of  the  type  being 

considered  is 

subjected  to  an  environment  where  a  uniformly  distributed  loading  of  the 

following  form  exists  in  a  transverse  direction  to  the  beam. 

' 

F(t)  =  F,sin7,t  +  ...  +  F  sin?  t  . 
1    1         n    n 

(3-D 

The  displacement  equation  for  this  beam  would  be  as  given  in  equation  (3-2). 

In  this  equation,  v.(x,t)  is 

n 
v(x,t)  =  £  vJ[(x,t) 

(3-2) 

i=l 

the  displacement  function  defined  by  equations  (2-1*0  when 

A'i2 

J     EI 

is  substituted  for  p.  If  any  of  these  P's  are  such  that  PL  is 

equal  to  or 

close  to  values  of  PL  in  appendix  C  list  1,  then  the  fixed  free  beam  will 

experience  large  amplitudes.  The  purpose  of  this  chapter  is  to  develop  a 

design  procedure  for  an  absorber  which,  when  connected  to  the  upper  end  of  the 

beam  will  decrease  these  amplitudes. 

Design  of  Single  Mass  Pendulum  Absorber 

Suppose  that  F(t)  in  equation  (3_l)  consists  of  a  single 

frequency  7  . 

Consider  a  one  mass  pendulum  which  has  a  natural  frequency  y. . 

If  7  is  the 

19 

20 


natural  frequency  of  the  pendulum,  then  it  follows  that  SL     must  be  as  defined 
in  equation  (3-3 )• 

',--%•  (3-3) 

71 

From  equations  (2-21),  it  is  apparent  that  D,  is  in  indeterminate  form. 
By  use  of  1'  Hospital's  Rule,  D,  becomes  as  defined  by  equation  (3"^b)  as 
7~»7   .  The  displacement  equation  of  the  beam  now  becomes  equations  (3"^)« 


F  sin/nt 

v(x,t)  =  — r 

P  EI 


_D1D2+D3j   ' 


(3-^a) 


_  2  coshftL  cosPL 
cospL  +  cosh31,  ^_         /,  r,\ 

sinPL  -  sinhBL  -  I  £7— THT"Jl  cosBL  -  coshpL 

\cosPL  +  coshpL/  \ 

'sinftL  +  sinhSlA  / 
J2       &XI1KA  "  axi'11'*-     VcosPL  +  coshPV  V 


D.  =  sinPx  -  sinh£3x-   ginpi,  +  smnoi,   ^^  _  ^^   _       (5^c) 


D3 


cosh3L  cosSx  +  cosBL  coshSx   ,  / ,  l  ,  \ 

L       cosBL  +  coshBL         J   "  * 

Notice  that  v(L,t)  is  equal  to  zero  in  equation  (3-i0« 

Since  cosBL  =f  coshBL  for  positive  values  of  PL,  D?  and  D,  will  always 
remain  bounded.  Therefore  for  the  absorbed  beam,  resonance  will  occur  at  y. 

only  when  the  denominator  in  equation  (3~^b)  is  equal  to  zero.  By  setting  the 
denominator  of  L\  equal  to  zero,  equation  (3"5)  is  formed. 

tanpL  =  tanh  PL  (3~5) 


* 

The  first  fifteen  roots  of  this  equation  are  contained  in  appendix  C, 

list  20 
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Equation  (3_5)  is  the  characteristic  equation  of  the  fixed -pinned  beam.1 
When  equation  ft-h)   is  evaluated  for 


0  = 


'pVi 


EI 


it  becomes  apparent  that  v(L,t)  =  0,  provided  equation  (3*5)  is  not  satisfied. 
Therefore,  if  equation  (3~5)  is  not  satisfied,  the  absorber  under  considera- 
tion effectively  "pins"  the  upper  end  of  the  beam  for  7=7,.     Resonance  will 
occur  at  7  =  71  if  and  only  if  7,  is  a  natural  frequency  of  the  fixed -pinned 
beam. 

The  pendulum  displacement,  v'(t),  as  given  by  equation  (2-21b)  is  in 
indeterminate  form  when  7  -   7\  •  By  use  of  1 'Hospital's  Rule,  equation  (2 -21b) 
is  reduced  to  equation  (3_6)  when  7  -  7-,,   and  tanPL  f   tanhPL. 


-pAF^siny  t  /  \ 

v/(t)  =  — =-* —  (  coshpL  sinpL  +  cospL  sinhpLJ  (3~6) 

1       Pp  EI  m    ^  / 

It  is  interesting  to  note  that  the  pendulum  displacement  is  inversely 
proportional  to  ra,.  Thus,  if  design  considerations  limit  the  pendulum 
displacement  to  small  amplitudes,  m.  must  be  larger  accordingly. 

It  should  be  noted  that  nothing  was  stated  about  resonances  or  lack  of 
resonances  of  the  absorbed  system  for  frequencies  other  than  71 .  The  complete 
spectrum  of  resonant  frequencies  for  the  absorbed  system  can  be  found  by 
solving  equation  (2-22). 


"TCin  N.  Tong,  Theory  of  Mechanical  Vibrations.  (New  York,  I963),  p.  258. 
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Design  of  Multiple  Mass  Pendulum  Absorber 

If  the  forcing  function  acting  on  the  beam  consists  of  more  than  one 
sinusoidal  component,  and  if  none  of  these  frequencies  are  such  that  equation 
(3~3)  is  satisfied,  then  it  is  possible  to  design  a  multiple  mass  pendulum 
absorber  which  will  effectively  pin  the  upper  end  of  the  beam.  Also,  the 
absorber  will  prevent  beam  resonance  from  occurring.  The  purpose  of  this  sec- 
tion is  to  show  how  such  an  absorber  can  be  designed. 

In  chapter  two,  it  is  shown  that  the  resonant  frequencies  of  a  canti- 
lever beam  with  a  multiple  pendulum  attached  are  roots  of  the  transcendental 
equation  (3~7)« 

0  =  D-l  +  D2   .  (3-7a) 

D  =  sinpL  -  sinhpL  -  (Wl  +  sinhgjA  f      _      \  <        (    ) 
1  \cospL  +  coshpL/  \  / 

-i  [EIp5(2  +  2  coshpL  cospL)] 

D2  = l .  (3-7c) 

n 

(7  m1   )  m.s..  J  f  cosPL+coshPL 

i=l 

D,  will  be  zero  if  and  only  if  equation  (3~5)  is  satisfied.  If  equation  (3~5) 

is  not  satisfied,  and  if  D?  is  forced  to  be  zero  at  the  forcing  function 

frequencies,  then  it  can  be  guaranteed  that  resonance  will  not  occur  at  these 

frequencies. 

2    V 

In  Appendix  B,  it  is  shown  that  7m    )  m.s.,  approaches  infinity  as  7. 

i=l 
approaches  7  (a  natural  frequency  of  the  uncoupled  pendulum).  Thus,  if  the  ab- 
sorber pendulum  has  as  its  uncoupled  natural  frequencies  the  frequencies  of  the 
forcing  function,  then  D„  will  be  zero.  Thus,  if  equation  (3~5)  is  not  satisfied 


23 


by  any  of  the  forcing  function  frequencies,  then  it  can  be  guaranteed  that 
resonance  will  not  occur  at  these  frequencies. 

In  order  to  design  a  multiple  pendulum  of  order  n  which  has  n  specified 
uncoupled  resonant  frequencies,  it  is  necessary  to  solve  equations  (3~8) 
simultaneously  for  the  pendulum  parameters . 

I  ~7±2W   -g(M)|  =0   *  (3-8) 

i  =  1,  • .  •  ,  n 

(M)and  (N)  are  defined  by  equations  (2-33)  and  (2-3*0  respectively.  An  n-mass 
pendulum  has  the  parameters  m  and  I , ,  i  =  1,  . ..,  n.  Thus  the  system  is 
underdetermined,  and  there  are  an  infinitude  of  suitable  pendulums  to  choose 
from.  Some  design  applications  may  place  restrictions  on  certain  masses  or 
certain  lengths.  Since  there  are  2n  pendulum  parameters  and  only  n  restraint 
equations,  considerable  freedom  is  left  to  the  designer  in  choosing  a  suitable 
pendulum. 

In  most  design  applications  only  a  few  of  the  frequencies  given  in 
equation  (3_l)  will  cause  the  unabsorbed  beam  to  have  large  amplitudes.  Thus, 
the  designer  may  feel  that  it  is  only  necessary  to  design  an  absorber  to  make 
the  coupled  system  non-resonant  at  these  frequencies.  Such  a  pendulum  design 
would  undoubtedly  be  easier  because  of  the  reduced  pendulum  order.  However, 
in  such  an  instance,  the  coupled  system  may  be  sensitive  to  frequencies  in 
equation  (3-1)  for  which  the  free  beam  was  unsensitive.  Thus,  it  is  best  to 
use  as  an  absorber  a  multiple  pendulum  which  has  all  of  the  frequencies  in 


It  should  be  noted  that  these  equations  are  nonlinear  in  m.  and  i.  . 
Hence,  a  numerical  solution  will  probably  be  necessary  when  n  is  large. 
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equation  (3~l)  as  uncoupled  resonant  frequencies. 

When  a  beam  is  acted  upon  by  a  forcing  function  as  given  in  (3~l),  and 
a  multiple  pendulum  absorber  (as  suggested  in  this  section)  is  attached  at  the 
upper  end  of  the  beam,  the  deflection  equation  is  somewhat  simpler  than  the 
deflection  equation  as  defined  in  equations  (2-k5).     The  simplified  displace- 
ment equations  become  equations' (3~9)» 

n 

v(x,t)  =  )  vt(x,t)  .  (3-9a) 

1*1 

{Vi(t)J   X  (u)i  Bln7it  *  (5"9b) 

i-1 


F  sin7,t  /BB 
vi(x,t)=-Sr^-f-P+B;  .  (3-9c) 


2  cos^L  cosh?  L 

Bl  =  1  "  cos^L  +  cosh^L  '  (5_9d) 

/sinP.L  +  sinhP  L\  /  \ 

B2  =  sinf^x  -  sinhP.x  -    (^^T^hP^ \cosV  "  coshp^j     .           (3"9e) 

/sinP.L  +  sinhp  L  \  /  \ 

h  =  Sin?iL  "  Sil^iL  "  VcosP^L  +  cosh^L  J  (,C0SpiL  "  COsh3iV     '            <5"9f > 


cosh?  L  cosp. x  +  cosp.L  coshP.x 
B4  =  cosPiL  +  coshp  L  •'"  1     '  (3"9g) 

{u)i  =  B5  B6  B?      .  (3-9h) 
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-F.pA 
P   ^i  EI 


coshP.L  sin£.L  +  cosP.L  sinh(3.L 

t>, i i i L_  fa-Qi) 

6~  cosP.L  +  coshf^L  u  >lW 


Ceo}. 
B^  =  __ 1 .  (5.9k) 

In  equation  (3~9k),  (cd)  is  the  eigenvector  of  the  uncoupled  pendulum  corre- 
sponding  to  the  eigenvalue  7.  .  cd  .  is  the  j   element  in  {co}   . 


PART  IV 
EXPERIMENTAL  PROCEDURES  AKD  RESULTS 


In  order  to  test  the  derived  equations  contained  in  previous  portions 
of  this  report,  various  experiments  were  devised  and  executed.  Since  experi- 
mental equipment  was  not  available  for  providing  a  uniform  sinusoidal  forcing 
function  as  considered  previously,  it  was  not  possible  to  test  the  displace- 
ment equations  as  such.  Instead,  a  point  forcing  function  was  used  in  the 
experiments  to  obtain  resonant  frequencies  of  various  cantilever  systems. 
Since  resonance  should  generally  occur  at  a  frequency  regardless  of  the 
distribution  of  the  forcing  function,  the  characteristic  equations  (2-15), 
(2-22),  and  (2-H^)  could  be  tested  by  using  a  point  forcing  function. 

A  description  of  the  parameters  involved  in  the  experiments  is  listed 
in  table  1.  In  experiments  B  and  C,  absorbers  were  designed  so  as  to  have 
uncoupled  resonant  frequencies  corresponding  to  the  first  resonant  frequency 
of  the  uncoupled  beam  (experiment  B)  and  the  first  two  resonant  frequencies  of 
the  uncoupled  beam  (experiment  C).  It  was  felt  that  this  would  be  a  more 
rigorous  test  of  absorber  performance  than  tests  in  which  absorbers  did  not 
possess  these  properties. 

Three  main  experiments  as  described  in  table  1.  were  conducted.  A 
schematic  of  the  experimental  apparatus  used  is  shown  in  figure  5.  Theoreti- 
cal resonant  frequencies  were  obtained  by  using  the  computer  programs  found  in 
appendix  D. 


26 


27 


TABLE  1 

EXPERIMENT  DESCRIPTIONS 

EXPERIMENT  A:  Unabsorbed  cantilever  beam 
BEAM  DIMENSIONS:   60"  X  l"  X  3/l6" 

BEAM  MATERIAL:  Reynolds  aluminum  No.  ,2021+  T351  QQA  268 
BEAM  PARAMETERS:  Afe  =  .188  in  ,  E  =  10. J  X  10  PSI 

p  =  .101  ~  ,  L  =  60"  ,  I  =  5.53  X  10"*  lnU 
in-5 

ABSORBER  DESCRIPTION:  No  absorber 

EXPERIMENT  B: 

BEAM  DIMENSIONS:  Same  as  experiment  A 

BEAM  MATERIAL:  Same  as  experiment  A 

BEAM  PARAMETERS :  Same  as  experiment  A 

ABSORBER  DESCRIPTION:  Single  pendulum,  I     =  3.5^",  n^  =  .101+  ibm 

EXPERIMENT  C: 

BEAM  DIMENSIONS:  Same  as  experiment  A 
BEAM  MATERIAL:  Same  as  experiment  A 
BEAM  PARAMETERS:  Same  as  experiment  A 

ABSORBER  DESCRIPTION:  Double  mass  pendulum,  I     =  2.00"  , 
/g  »  I.63",  m  =  .028  lbm,  m2  -  .251+  Ibm. 
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The  theoretical  and  experimental  resonant  frequencies  for  the  various  experi- 
ments are  listed  in  table  2. 

The  experimental  resonant  frequencies  were  found  by  starting  veil  below 
the  first  resonant  point,  and  slowly  increasing  the  frequency  until  beam  dis- 
placements increased  substantially  over  a  small  frequency  range.  The  fre- 
quency was  then  increased  until  displacements  decreased.  The  frequency  range 
of  large  displacements  was  then  retraced  at  a  slower  rate  in  order  to  determine 
the  frequency  at  which  displacements  were  largest.  Since  the  frequency  range 
contained  large  displacements,  the  amplitude  of  the  forcing  function  was 
reduced  in  order  to  preserve  the  linearity  of  the  system.  An  electronic 
counter  was  then  used  to  determine  this  frequency. 

The  lover  resonant  frequencies  were  very  easily  detected.  At  higher 
frequencies,  however,  resonant  detection  was  more  difficult  due  to  smaller 
displacements.  For  this  reason,  only  the  first  several  resonant  frequencies 
were  found  experimentally. 

In  general,  the  experimental  resonant  frequencies  were  about  10jo  higher 
than  the  theoretical  resonant  frequencies.  It  was  felt  that  the  main  reason 
for  this  discrepancy  was  a  slight  error  in  the  value  of  E  used  in  the  analyti- 
cal work.  If  a  larger  value  of  E  had  been  used  in  the  analytical  work,  then 
the  results  in  experiment  A  would  have  been  closer  in  agreement.  The  dis- 
crepancies in  the  results  of  experiments  B  and  C.  were  probably  caused  by 
slight  errors  in  the  values  of  E  and  the  pendulum  parameters  used  in  the 
analytical  solution.  It  vas  felt  that  the  agreement  of  results  in  table  2  was 
sufficient  to  validate  the  characteristic  equations  of  chapter  two. 

The  first  several  resonant  mode  shapes  for  the  cantilever  systems  con- 
sidered were  determined  by  using  the  computer  programs  found  in  appendix  D. 
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Normalized  mode  shapes  for  experiments  A,  B,  and  C  are  shown  in  figures  6,  "J, 
and  8  respectively. 
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PART  V 
SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  shown  that  for  a  beam  which  is  forced  transversely  by  a 
uniform  forcing  function  F  =  F-^in?  t  +  ...  +  F  sin?  t,  a  properly  designed 
multiple  pendulum  absorber  attached  at  the  unfixed  end  will  prevent  resonance 
from  occurring  under  most  conditions.  If  the  multiple  pendulum  has  as  its 
uncoupled  resonant  frequencies  the  frequencies  of  the  forcing  function,  then 
lack  of  resonance  can  be  guaranteed  provided  none  of  the  7.,'s  correspond  to 
natural  frequencies  of  the  fixed-pinned  beam.  Under  these  conditions,  it  can 
also  be  guaranteed  that  the  unfixed  end  of  the  beam  will  be  effectively  pinned. 

Since  resonant  frequencies  are,  in  general,  properties  of  the  physical 
system  and  not  properties  of  the  forcing  function  distribution,  the  lack  of 
resonance  can,  in  general,  be  guaranteed  even  if  the  forcing  function  is  not 
uniform  along  the  length  of  the  beam. 

In  design  instances  where  the  forcing  function  contains  a  large  number  of 
frequencies,  the  design  of  a  multiple  pendulum  of  this  high  order  may  not  be 
feasible.  An  alternative  approach  would  be  to  design  a  multiple  pendulum  (of 
somewhat  lower  order)  which  had  as  its  uncoupled  natural  frequencies,  only  the 
forcing  function  frequencies  with  large  coefficients.  Lack  of  resonance  would 
then  be  guaranteed  at  these  frequencies.  The  resonant  frequencies  of  the 
coupled  system  could  then  be  obtained  by  using  the  computer  programs  found  in 
appendix  D.  If  none  of  these  frequencies  corresponded  to  frequencies  of  the 
forcing  function,  then  lack  of  resonance  could  be  guaranteed  for  the  coupled 
system. 
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APPENDIX  A 

SUPERPOSITION  PROOF 

It  vas  shown  that  for  a  uniformly  distributed  forcing  function  with  a 

single  frequency,  the  following  solution  exists  for  the  displacement  equations. 

v(x,t)  =  r(x)sin7t  . 

v'(t)  =  {ujsinrt  . 

It  will  now  be  proven  that  if  the  forcing  function  F  =  F^inr-jt  +  ...  + 

F  sin?  t  acts  on  the  beam,  then  the  resultant  displacement  equations  will  be  as 
r    r 

shown  below. 

v(x,t)  =  \  rjLsinrit  .                    (A-l) 

i=l 

r 

v'(t)  -  V  {u)isin7it  .                     (A -2) 

1-1 

In  the  above  equations,  r.sin7.t  and  {u^sin  7±t   are  the  resultant  solutions 

when  the  forcing  function  is  F  =  F^inr^. 

The  governing  differential  equation  and  the  boundary -conditions  for  a 

multiple  frequency  forcing  function  are  given  below. 

o\(x.t)   ^  a2v(x,t)      y  Fi5lnv                       (A.5a) 
^     EI     ^2     "  4    EI      ' 

v(0,t)  =  0  .                                            (A-5b)       j 

dv(x,t) 

=   0  .                                        (A -3c) 

x=0 

dx 
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d2v(x,t) 

at2 

=  0  .                                    (A -3d) 

x=0 

♦ 

EI  jSji^SJ 

A   A   a2v-  (t) 

ox     x=L   i=l   3=1 

It  will  now  be  shown  that  equations  (A -3)  are  satisfied  when  v(x,t)  and  v'(t) 

are  as  defined  in  equations  (A-l)  and  (A-2). 

By  substituting  equation  (A-l)  into  equation  (A -3a),  equation  (A-ka)   is 

obtained . 

o\(x,t)   pAfe   ^  *\(x,t) 


F,sin7,t 


i=l   dx  1«1   dt       i=l 

+ v» 
This  equation  is  obviously  true  because  it  is  true  for  each  i   component.  By 

substituting  equation  (A-l)  into  equations  (A -3b),  (A -3c),  and  (A -3d),  equa- 
tions (A -Ub),  (A-Uc),  and  (A-lid)  are  obtained. 


r 
£  vi(0,t)  =  0  . 


1=1 
r 


I  &».*> 


i=l 


=  0  . 


(A -Jib) 


(A -lie) 


x=D 


*-  a2 


Itt^)  U  =  0  • 


(A -lid) 


i=l 


Note:   .v  (t)  is  the  displacement  of  m  due  to  the  loading 


F .sin?  t  . 
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As  In  the  case  of  equation  (A-Ua),  these  equations  are  valid  because  they  are 
valid  for  each  i   component. 

By  substituting  equations  (A-l)  and  (A -2)  into  equation  (A-Je)  the  last 
boundary  condition  becomes  equation  (A -lie). 

i-l     °X       i=l  q-1 

As  in  the  previous  cases,  this  equation  is  valid  because  it  is  valid  for  each 

1th  component.  Thus  the  resultant  displacement  equations  due  to  a  loading 

r 
function  of  the  form  F  =  L    F  sin  y  A,   are  as  shown  in  equations  (A-l)  and 

i=l 
(A-2). 
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APPENDIX  B 
Singularity  Limit  Proof 

The  purpose  of  this  appendix  is  to  prove  that  as  J   approaches  7  (where 

7  is  a  natural  frequency  of  the  uncoupled  pendulum),  then  {s  , }  approaches 
s  l-i- 

«{a>}.  {«>}  is  the  normalized  eigenvector  of  the  uncoupled  pendulum  correspond- 
ing to  7  .  This  proof  is  divided  into  3  parts. 
s 


Part  1 


The  purpose  of  this  part  is  to  show  that  the  two  multiple  pendulums  in 
the  following  figure  have  no  common  resonant  frequencies. 


\\  N 


■\\N 


6 


6 


-pendulums  are  identical 
below  this  line 


6 


Figure  9* 

Similar  Multiple'  Pendulums 
These  two  pendulums  are  identical  except  for  one  added  mass  and  link  as  shown. 

When  pendulums  1  and  2  are  two  mass  and  one  mass  pendulums  respectively, 

the  resultant  characteristic  equations  are  as  given  in  equations  (A-?)  and 


kl 

(A-6).                                                             j 

p    gnu 

ONE  MASS  PENDULUM:   0  -  -7  m0  +  —■     .                              (A -5) 

TWO  MASS  PENDULUM: 

0  = 

2      mi  2   m2  2«          m2  2 

. 

(A-6) 

-g  m22            2     g  m22 

£                               ~7   m2  +  '     I 
n               2     ^2 

• 

Notice  that  the  lover  right  element  in  equation  (A-6)  is  the  right  side  of 

equation  (A -5).   If  these  two  pendulums  have  a  common  natural  frequency,  then 

the  common  frequency  must  be  as  defined  by  equation  (A-5).   If  this  is  so, 

/gm2\2 
then  by  equation  (A-6),  (  .     )    must  be  zero.   Since  this  is  impossible,  the 

two  pendulums  do  not  have  a  common  natural  frequency. 

When  pendulums  1  and  2  are  3  mass  and  2  mass  pendulums  respectively, 

then  the  resultant  characteristic  equations  are  as  shown  in  equations  (A -7) 

and  (A -8). 

2      m2  5    "S  5            g  "S  3 

(A -7) 

3                    ^3 

• 

k2 


0  = 


2     ^1  3 

'7ml  +  ~L 


gm, 


2  5 


gm 


2_1 


gm. 


2J. 


>    7  m  +  _2JL 


gm 


5_i 


gm 


±1 


gm 


2JL 


3 

gm 


>  -y  ^  + 


u. 


A-8) 


Notice  that  the  elements  in  equation  (A -7)  are  identically  equal  to  the  2X2 
array  of  elements  in  the  lower  right  corner  of  equation  (A -8).  By  denoting 
the  characteristic  determinants  by  D  ,  Dp,  and  D  for  a  3  mass,  2  mass,  and 
1  mass  pendulum,  respectively;  equations  (A -7)  and  (A-8)  become  Dp  =  D  =  0 
when  equations  (A-7)  and  (A-8)  are  satisfied  simultaneously.  By  expanding 
equation  (A-8)  by  minors  of  the  first  column,  equation  (A -9)  is  obtained. 

m„  ,\2 


0 


•7m1  +  g-^+g-f-^D2 


2_1 


D„ 


(A -9) 


If  equations  (A-7)  and  (A-8)  are  satisfied  by  the  same  value  of  7,  then 

D  =  Dp  =  0.  From  equation  (A-9),  D  =  0  for  this  value  of  7.   I  previously 

proved  that  D,  and  Dp  are  not  both  zero  for  any  common  value  of  7.  Hence, 

the  3 -mass  and  the  2-mass  pendulums  under  consideration  have  no  common  natural 

frequencies. 

It  will  now  be  proven  that  if  a  K-l  mass  pendulum  and  a  K  mass  pendulum 
have  no  common  natural  frequencies,  then  a  K  +  1  mass  pendulum  and  a  K  mass 
pendulum  will  have  no  common  natural  frequencies.   It  must  be  kept  in  mind 
that  the  pendulums  under  consideration  are  identical  except  for  the  addition 
of  masses  and  lengths  at  the  upper  ends  of  the  pendulums.  The  preceeding 


hi 


statement  vill  be  proven  by  contradiction. 

Assume  that  the  K-l  mass  pendulum  and  the  K  mass  pendulum  have  no  common 
natural  frequencies.  Assume  also  that  the  K+l  mass  pendulum  and  the  K  mass 
pendulum  do  have  at  least  one  common  natural  frequency.  The  characteristic 
equation  for  the  K+l  mass  pendulum  is  given  in  equation  (A -3.0). 


0  = 


2        L  ml  k+l       »": 

'7  rai +  ii    *  ~t 


2  k+l 


gm, 


2  k+l 


0 

t 


gm2  k+l 


th 


0  ,...,   0 


k   order  characteristic 
determinant 


k-l   order 

characteristic 

determinant 


(A -10) 


By   expanding  equation  (A-10)  by  minors  of  the  first  row,  equation  (A-ll)  is 
obtained . 


f  s   n  k+i  A  m2  k+i\  _  '  (m2  k+iN 


Dk  "  \~T~J  Dk-i    • 


(A -11) 


If  D.  n  and  D,  are  both  zero  for  some  value  of  7,  then,  by  equation  (A-10), 

D  n  must  be  zero  for  the  same  value  of  -/.     Since  it  was  assumed  that  D,  ,  and 
k~l  K~X 

D,  are  not  both  zero  for  any  value  of  y,   D,  -  and  D,  are  not  both  zero  for  any 

value  of  7. 

It  was  shown  previously  that  a  3~*rass  pendulum  and  a  2  mass  pendulum  of 

the  types  being  considered  have  no  common  natural  frequencies.  Hence,  by 

using  the  previous  argument  and  mathematical  induction,  the  n  mass  pendulum 

and  the  n-1  mass  pendulum  being  considered  have  no  common  natural  frequencies. 


hh 


Part  2 


The  purpose  of  this  part  is  to  show  that  as  7  approaches  one  of  the 

values  of  7  then  s,,  approaches  °°  .  s,,  is  the  first  element  in  the  column 
s      11  11 

vector  {s.,},  from  equation  (2-38), 

(Adj  r  ) 
(S)  =  (R)  1   -  rrM-     •  (A-12) 


The  element  s,,  now  becomes: 

(-1)2  D  , 

"n  =■ iT^  •  (A-15) 

11      Dn 


*  From  part  1,  D    does  not  approach  zero  as  D  approaches  zero.  Since  D 
approaches  zero  as  7   approaches  7  ,   s,,  approaches  infinity  as  7   approaches  7  « 

Part  3 

The  purpose  of  this  part  is  to  show  that  {s.,}  approaches  °°(u>}  as  7 

approaches  7   .  From  equation  (A-12), 
s 

{Adj  r.J 
{s   j  .  l_il_  (A.lU) 

lsilJ      R 


I  shall  now  solve  for  {co}  ,  assuming  that  the  first  element  in  {to}  is 


unity. 

(D,  =  1  .  (A -15a) 

(R)(a)}   =  {0}  .  (A-15b) 

s 

Expanding  equations  (A-15)  yields  equations  (A-l6)  . 


^5 


rl2a)2+  ...  +rlncon  =  -rn 


r22co2+  ...  +r2n%  =  -r21 


(A-16) 


n2  2        nnn    nl 
Solving  the  lower  n-1  equations  of  (A-l6)  by  Cramer's  rule  results  in  (A-17). 

'in 


r22   '   "rll  '        '>  r' 


r23   '   "r21  '  "*'  r2n 


2n   '     nl 


adj   r 


11 


nn 


(A-17) 


The  i-1  column  in  the  numerator  in  equation  (A-17)  is  the  right  side  of 
(A-l6). 

Noting  that  (R)  is  symmetric,  and  that  the  interchange  of  any  two 
columns  in  a  determinant  results  in  a  sign  change;  it  becomes  apparent  that 
equation  (A-17)  is  equivalent  to  equation  (A-l8). 


(-1)1*1  Mil    Ad'  ril 

to  =  ——————  c    . 

i    Adj  r„      Adj  rn 


(A-18) 


11 


11 


By  comparing  equations  (A-lU)  and  equation  (A-l8);  it  becomes  apparent 

that  as  7   approaches  7  ,   and  (R)  approaches  singularity,  then  (s  }  becomes 

s  ij_ 

as  defined  in  equation  (A-19). 


lim  {s.,}  =  ~  (co) 


(A-19) 


n 


It  is  interesting  to  note  that  in  some  instances  L    a> .    =   0,  but  in  no 


instance  is  equation  (A-20)  satisfied. 


i=l 


k6 

n 

^  ©jBij  =  0  .  (A -20) 

i=l 

If  equation  (A-20)  were  true,  the  momentum  of  the  pendulum  in  the  horizontal 

direction  would  necessarily  be  a  constant.  This,  in  turn  would  mean  that  l 

would  have  to  be  vertical  at  all  times.  If  /..  were  vertical  at  all  times,  /? 

would  have  to  be  vertical  at  all  times,  etc.  Thus,  the  trivial  case 

(co)  =  (O),  would  exist. 

n  n 

Since   )  m.u>  is  never  zero,   lim   7m    ;  m.a.-  «  «  .         (A-2l) 
i-1  7~*  7s  i»l 


47 


APPENDIX  C 
ROOTS  OF  TRANSCENDENTAL  EQUATIONS* 


Root  LIST  1  LIST  2 

No.  


COSHpL  COS3L  +1  =  0  TAHpL  -  TANHpL  =  0 

1  1.87  3-93 

2  4.69  7.07 

3  7.85  10.21 

4  10.99  13.35 

5  14.14  16A9 

6  17.27  19.64 

7  20.42  22.78 

8  23.56  25.92 

9  26.70  29.06 

10  29.84  32.21 

11  32.98  35.35 

12  56.12  38.49 

13  39.26  41.63 

14  42.41  44.77 

15  45.55  47.91 

*  jt 
This  is  only  a  partial  list.     For  n  >  15,   Root  list  1  ^nit  -  §■ 


2 


Root  list  2  ~  nit  +  r     • 
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APPENDIX  D 
COMPUTER  PROGRAMS 


inniiuniinii/iiifiiummmminijiiii/imnjiitj 


£::::::::::::*:n^n:**n********^^ 

C   SUBPROGRAMS  NEEDED  ARE...  T1AC,  vc  v  v^  yi  IPRIKT1 

C   SUBROUTINE  ROOT  (  F  ,  DEL  I  ,  DDEL ,  DELX ,  XMAX,  I  L  AST ,  XS,  Xu  ,  YS  ,  YL  ,  I  PR  n. .  ) 

C       FUNCTION    FUN  CXI.        .-•        ,ftM    C|lfc(fv.     fnBRF„nNI)ING    T0    THE    OPTION    DESIRED 


101  F0RMAT(AE14.5) 
DELI=.01 

X=.5 

■  CDEL=.l 
D£LX=.01 
XMAX=100. 
ILAST=1 
IPRINT=0  . 

CALL0ROOUX, DELI, DDEL, DELX, XMAX, ILAST, XS, XL, YS, YL , I  PRINT ) 

X=XL+3.05 
WRITE{3,101)XS,YS,XL,YL 

100  CONTINUE 
STOP 

unuminninnnuininnminniunnniininnumimin 

C   FOR  A  CANTILEVER  6EAH  WITH  A  SINGLE  MASS  .PENDULUM  ATTACHED  A,  THE 

C   UPPER  END. 

C   SUBPROGRAMS  NEEDEO  ARE....... 

C   FUNCTION (F)   OPTION  3 


101  FCRMAT(AE14.5) 


^9 


DELI-.l 

F  =  .2 

DDEL=.i 

DELX=.01 

XMAX=10000. 

ILAST*1 

IPRINT=0 

001001=1,20 

CALL    ROOT(F,DFLI,DDEL,DELX,XMAX,ILAST,XS,XL,YS,YL,IPRlNT} 
F  =  XL 

WRITE(3,101)XS,YS,XL,YL 
100    CONTINUE 
STOP 
END 
C//// /////////////// ///////////////////////////////////////////////// 


C/////V//// //////////////////////////////////////////////////// ////// 

C       THIS    IS    THE    MAIN    PROGRAM    FOR    CALCULATING    THE    DISPLACEMENT    EQUATION 
C   FOR  A  VERTICAL  CANTILEVER  BEAM  WITH  A  SINGLE  MASS  PENDULUM  ABSORBER. 
C  *  *  *  v  *  *  *  *  *  *  *  *  *  ********  *  *  if  *  *  *  *  *  *  *  *  * *  *  $  $  *  *  $  *  *  *  *  $  *  $  *  *  *  *  *  *  *  *  *  £  *  *  *  *  *  *  *  *  *  *  *  $ * 

£*******  *  *  *  *  *  *  *  $  if  $  *  X:  *  *  *  *  *  *  *  *  * *  *  *  *  $  *  *  *  #  *  *  *  *  *  *  *  *  *  *  *  *  *  $  *  *  *  V  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

C       BEAM    DISPLACEMENTS    FOR    THE    UNABSORBED    BEAM    ARE    ALSO    CALCULATED    FOR 
C       THE    SAME    FREQUENCIES    AND    VALUES    OF    X. 

Q  *  *  >;< *  *  *  $  *  *  *  * >;c  *  * *  $  *  *  •-}  *  #  *  *  *  *  * *  *  *  *  * *  *  *  *  *  *  *  *  * *  *  * *  *  *  *  *  *  *  *  *  *  * *  *  *  *  *  $  *  *  $  *  $  *  *  $  *  * 
C*+<t*  *  **************  *************************************************** 

C       THE    BEAM    FOR    WHICH    THE    PARAMETER    ARE    ENTERED    IS    THE    BEAM    IN 
C      EXPERIMENTS    A,B,AND    C. 

C  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  v  ft  *  *  *  ******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 
C***JJ  =)  ***************************************************************** 

C      THE    BEAM    LOADING    IS    1.0    POUND    PER    INCH. 

£*****  *  *  *  *  *  *  *  *  *  *  $  <:  *  X:  *  *  *  *  *********  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  ******  *  *  *  *  *  *  $  *  * 
C***$*  ***************************************************************** 

C       BEAM    DISPLACEMENTS    ARE    CALCULATED    FOR    X=0, 6, 12, 18, . . . , 60    INCHES. 

Q****  *  *********************************************** *******i;;*** ******* 
C****<  ***************************************************************** 

C      THE    FREQUENCY    RANCE    IS    1.0    RAD. /SEC.    TO    300.0    RAD. /SEC.       AS    THE 
C       PRCGRAM    STANDS    . 

C  *  *  *  *  *  *  *  *  *  *  *  *  **************  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  if  if  *  *  if  if  if  if  if  if 

c***** ***************************************************************** 

C       SUBPROGRAMS    NEEDED    ARE 

C   FUNCTION  SUBPROGRAM  VF8(X,B,AL) 

C       FUNCTION    SUBPROGRAM    VSM < X , B, AL , E* A  I , AM, D, G ) 

C  *  <:  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  if  *  *  * 

c**** ************** **************************************************** 
c*** ***************************************************** ************** 

C    A.^PSNDULUM  MASS  IN  LBM. 

C    C=PENDULUM  LENGTH  IN  INCHES. 

C    AL=EEAM  LENGTH  IN  INCHES. 

C    E=MODULUS  OF  ELASTICITY  IN  PSI. 

C         A3=CR0SS    SECTIONAL    AREA    OF    BEAM. 


_ 
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C  AG=ACCELERATION    CF    GRAVITY    (386    IPS). 

C  ROBEAM    MATERIAL    DENSITY    IN    LBM/IN**3. 

C  6=BETA. 

C  FR=FREQUENCY    (CPS). 

C  G=GAMMA       (RAO/SEC). 

C  VF(J)=FREE    BEAM    DISPLACEMENT. 

C  VP( J)=ABS0RBED    BEAM    DISPLACEMENT. 

C  V?P= PENDULUM    MASS    DISPLACEMENT. 

Q  $  $  if  £  if  $  $  i-.  if  it  if  if  if  if  #  if  $$«$#$$#  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  <:  *  *  ;: 
f  «  if  *  *  *  «  *  *  *  *  «  «  #  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  #  #  *  *  *  *  *  *  *  *  *  ###**#*#####  *  *  *  ♦  *  *  *  *  *  *  *  *  *  *  *  *  *  * 
Q  *  #  V  «  <  $  S)s  «  *  if  if  $  if  «  *  *  $  *  *  *  * ){:  *  *  <<  «  *  #  * «  if  $  #  *  *  if  if  if  if  $  *  ##*>•,(##  <<jCo;c>:<  #  *  #  *  #  $  $  #  i(c*  *  4  *  #  *  if  if  if  if 

DIMENSION    VF(20),VP(205 
10 i    F0RMAT(1H0,F7.3,F7.3,F6.5,11E10.3) 

102  FORMAT (1H    ,  F7. 3, F7. 3, F6. 5 , 1 1E10. 3) 

103  FORMAT (LH    ,E10.3) 
AM=.i04 

D=3.54 

AL=60. 

E=10.3E6 

AI=5.53E-4 

AB=.188 

RO=. 101 

AG=386. 

0011=1,3000 
^      G=FL0AT(I) 

FR  =  G/2. 

B=( {RC*AB*G*G)/(E*AI*386. ) )**.25 

X=-6. 

GS=G*G 

D02J=1,U 

X=  X  +  6  . 

VFi  J)=VFBU,B,AL) 

VPf JJ=VSM(X,3,ALtE,AI, AM,OtG) 
2    CONTINUE 

WR 1 TE { 3 , 1QL ) G, FR, B, C VF  (J ) , J= 1 » 1 1 ) 

WRITE(3,102JG,FRtB, (VP( J) , J=l,ll) 

VPP=VP(ll)/( i.-GS*D/AG) 
1    v:RITE(3,103)VPP 

STOP 

END 
C/  / 7  //  / 7 7 7 '/ V I 7 '/// 7 7 1 7 '// 7/ 7 '// '/// 7 7/ '// 7 '// 7/ 7/ 7 7 '// 7 7 '// 7 7 7 7/ 7/ 7/ 7/ 7/ V 7/ 7  / 
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Q*****************************************^*****^ 

c 

C*  = 

C*> 


DEFLECTIONS    OF    UNABSORBED  -BEAM    ARE    ALSO    CALCULATED.  ...... 

**************************  *********#**************************vik:v'':sr'!f: 

**** * *********** *#«************#***************************** v*7**v** 


c***< 

c 
c 

c 
c 
c 

t  ^BEAmIaRAMETERS  ARE  THOSE  USED  IN  EXPERIMENTS  A, B.AND  C. 

t**^***********************************^ 

^^<ty:^^<^*«,:*******************************':=*********::t****):t*,:!:;!)^*':;^^ 
C   PENDULUM  PARAMETERS  ARE  THOSE  USED  IN  EXPERIMENT  C.  

j£************************************^ 

C  SUBPROGRAMS  NEEDED  ARE 

C  FUNCTION  VFS(X,3,AL) 

C  SUBROUTINE  VMP ( X , B , AL , E , AI , G , Dl , VNP ) 

C  SUBROUTINE  MI NV ( A, N , 0, L , M) 

£******^*^**^ 

C  AM(I)=MASS  OF  THE  ITH.  PENDULUM  MASS   (LBM.) 

C  CL(I  ) -LENGTH  OF  LINK* I  ) 

C  K=CRDER  OF  THE  PENDULUM 

C  VF(J)=FREE  BEAM  DISPLACEMENT  Al  X  =  6*J  .(IN) 

C  VAB(J)=ABSCRBED  BEAM  SISPLACEMENT  AT  X=6*J   (IN) 

C  VP?(L)=DISPLACEM6NT  OF  LTH.  PENDULUM  MASS  IN  INCHcS 

C    AL=3EAM  LENGTH  IN  INCHES. 

C    E=MODULUS  OF  ELASTICITY  IN  PSI. 

C    AB=CROSS  SECTIONAL  AREA  OF  BEAM. 

C  AG=ACCELERATION  OF  GRAVITY  (336  IPS). 

C  '  RC=BEAM  MATERIAL  DENSITY  IN  LBM/IN**3. 

C    B=BETA. 

C    FR=FREQUENCY  (CPS). 

G    GsGAKMA   ( RAD/SEC ) .  .„...*.  v-  •  - &*i 

c***,**«**********************^^ 

r************************************^^ 

c******************************^***<:**********************^***^**vvvv 

DIMENSION  VF(20),VPP( 10) ,VAB(20) 
COMMON  AM(IO) ,CL(L0),AR(100),K 

102  FORMATdH  ,  F7.  3  ,  F7  .  3  ,  F6.  5  ,  1 1EL0.3  ) 

103  FORMATdH  ,  E1C.3) 
iGA  FORMAT (10c 12. 3) 
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105  FORMAT (15) 

1C9  FORMAT ( I 7H-START  OF  NEW  RUN) 

K=2 

AM! I)  =  .02 76 

AM(2)=.2535 

CU1)=2.0 

CL!2)=1.635 

AL=60. 

E=1G.3E6 

AI  =  5.53E--4 

A3=.188 

RO=.101 

AG=386. 

0011=1,300 

G=FL0AT(U/4. 

FR=G/2./3. 14159 

B={  (R0*AB*G<CG)/(E*AI*386.  )  )**.25 

X=-6- 

GS  =  G*G' 

D02J=l7il 

X=X+6. 

VF(J)=VFB{X,B, AL) 

GALL  VMP\X,B,AL,E,AI,G,D1,VNP) 

VAB(J)=VNP 

2  CONTINUE 

WR I TE ( 3, 102 ) G, FR, B, ( VF ( J ) , J=l» 11 ) 
WRITE (3,102>G,FR,3, (VAB(J) ,J=1 til) 
D03L=1,K 

3  VPP(L)=(AR(L) )*VAB(11)*(AAM(1) )*AG/(CL(1) ) 
WRITE (3,104) (VPP(L) ,  L=1,K) 

1  CONTINUE 
STOP 
END 

ummmiuunuiiimmumumu/titmumiuuiimmiii 


SUBPROGRAMS 


zuiiiuntniiniiiiimiuiiui/nmi/niimuini/mimmn/i 

FUNCTION  FON(X) 

£#**«$  $$  $«***£****************************  *##**#***#***#♦******#******#' 

C   OPTION  1 

r****4 ********«*##***********************************************#***** 

FUN=(COSH(X) )*(COS(X)  )  +  l. 
RETURN 
END 
C/ V/  /  /  // I 1 II 7 7 7 7 7 7 7 '/  // 7 7 7 7 7// 7 7/ '// 7/ 7 7 7 7 7 7 7 7 7 7/ 7 '// 77 7/ 7 7 7 7 7/ 7/ 7/ 7 I II 7 
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c///////////////////////////////////////////////////////////////////7 

FUNCTION  FUN(X) 

Q*****  **********  ******************************************************* 

C   OPTION  2  .  

£«***  4  *************  **************************************************** 

FUN=TAN(X)-TANH(X) 

RETURN 

z/  / 1 1  m  1 1 1 1 1 1 1  u  /  n  j  i  u  u  i  m  m  i  /  u  n  u  i  n  n  1 1  u  i  /  u  n  i  n  m  i  /  n  //  / 


ujjiiiniiJi/JUiiiiiiii/ii/ni/ninjinfin/ff/n/n/u/mnnin 

FUNCTION    FUN(F)  „.,.....,. 

£*************************:******************^^ 

c*******'*********************^ 

G=F*2.*3. 14159 

GS=C*G 

R0=.101 

AL=60. 

E=10.3E6 

AI=5.53F-4 

AB=.188 

AG=386. 

B4=RO*AB,*G*G/(E*Al*386.) 

B=B4**.25 

Bl=B*AL 

AM=.10A 

0=3.54 

011=1. /GS-D/AG 

D1  =  -AM/DU 
D2=E*AI*B*B*B/{COS(BL)+COSH(BL) ) 

02=02*386. 
C2i=2.*2.*C0SH{BL)*C0S{BL) 

C22=D2*C2I 
C23=(S:n<BL)+SINH<BL))*(CGS(BL>-COSHCBL)) 

G24  =  C23/(C0S(eD+C0SK(BL)  ) 

C2  5=SIN(EL)-SINH(BL)-C24 

C2=D1*C25+C22 

FUN=C2 
RETURN 

C// 7 7 '//// 7 7 7/ V  // V V/ 7 7 7/ 7 7/ 7/ 7/ 7 7/ /// V/ 7/ 7/ /// 7 7 7 7/ 7 7 7 7 '// 7/ 7/ 7/ 7 7 7 '/// 7 


C/  /  /  /  /  //  /  /  /  / 7 7 7 7 7 7/ 7 7/ 7 7 7 7 7 7 7 7 7 7 7/ 7/ 7/ 7 7 7/ 77 7/ 7/ 7 7 7/ 7/ 7 7/ 7/ 7 7 7 7 7 7/ 7 7/ 

SUBROUTINE  ROOT ( X , DEL  I , DDEL , DELX , XMAX , I  LAST , XS , XL » YS7 VL , I  PR  I  NT ) 
C   SUBROUTINE  ROOT  WAS  WRITTEN  BY  OR.  H.  S.  WALKER,  MECHANICAL  ENGINES: 
C   DEPARTMENT,  KANSAS  STATE  UNIVERSITY 
C      THIS  SUBROUTINE  CALLS  FUN(X)  *********************************** 


:?4 


c 
c 
c 

u 
C 

c 
c 


10 


98 

99 

100 


I  PRIM 

ROOTS 

DELI 

CDEL 

OELX 

XMAX 

ILAST 


1=1    PRINTS    XS,XL,YS,YL.       IPRINT=0    ODES    NOT    P 

OF    THE    TRANSCENDENTAL    EQUATION    FUN(X)    ***** 

IS    THE    INITIAL    INCREMENT    ******************* 

FACTOR    FOR    DIVIDING    THE    INTERVAL    **** 

ACCURACY    ON    THE    SOLUTION    *MMMMMMMMMM 

MAXIMUM    THAT    THE    ARGUMENT    X    MAY    BECOME    * 
NUMBER    OF    ROOTS    DESIRED    *************** 


IS 

IS 
IS 
IS 


THE 
THE 
THE 
THE 


-,-  .1  V  *  V  >r  V  -i-  f  ' 

if  if  if  if  if  if  if  if  if  ■ 


FORMAT (AE20. 8) 

IUI  =  0 
X  =  X-0 
DEL  = 
X  =  X 
Y=FUN 
IF(Y) 
ZAZ=- 
GO  TO 
ZAZ  =  + 
XS    = 


110 
120 
130 

140 


YS  = 
IF(X. 
X  =  X 
Y  =  FUN 
I  F  (  Z  A 
XL  = 


ELI 

DELI 

+    DEL 
(X) 

98,500,99 
1.0 

100 
i.O 
X 

V 

* 

CT.XHAX)       GO    TO    510 

+    DEL 

:x) 

Z*Y)     110,500,100 


L    =    Y 

=  XS  - 
F(XL-XS 
=FUN(X) 
EL  =  DE 
F(ZAZ*Y 
L  -  X 
L  =  Y 

=  X  - 
=FUN(X) 


YS*(  (XL-XS)/(YL-YS)  ) 
-OELX)    500,130,130 

L*DDEL 

)     140,500,100 


150 


500 


900 
502 

510 


DEL 

IF(ZAZ*Y)     140,500,150 

XS    =    X 

VS    =    Y 

SO    TO    12 

I U I  =  I U I  + 

IF  I I  PR  IN 

RITE(3, 

ONTINOE 

F  i I  LAST 

ONTINUE 
*    *    * 


0 

i 

L 

T.EO.O)    GO    TO    900 
4)   XS,XL,YS,YL 

-IUI)  510,510,10 
*********************** 


V   5?   * 


RETURN 
END 
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c////  /  / in i tin  mu  mi i n iiu/uu imiimmmmi '////// ////// // 

FUNCTION  VF3(X,BfAL) 

C      DEFLECTION  FOR  BEAM  WITH  NO  PENOOLUUM  . 

C  *  #  *  *  >>  «  *  *  *  $ *  *  #  *  *  ******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

(;**#>>>}  ****  *********  4  *************************  ************ ************* 

AI=5.5E-4 

E=10.3E6 

Bi*8*AL 

P  *'  —  £;  &  Y 

C  1  =  1  COSH  {  3D  )  *  (  S  I N  (  BL)  )  *  I  COS  (  BL  )  )  *  {  S  I  NH  {  BL  )  ) 

C2=Cl/(2.+2-*(C0S(8L) )*(COSH(BL> ) ) 

C3=COS(BX)-COSH(BX) 

C4=C3*(SIN(8L)+SINH{BL) 5 / ( COS { BL ) +COSHJ BL ) ) 

C5=SIN(BX)-SINH(BX)-C4 

C6=(C0SH(BL) )*(COS(BX))+(COS(BL))*(COSH(BX) ) 

C7=C6/{C0S(BL)+C0SH{BL) ) 

C8-C7-1. 

VFB=(C2*C5+C8)/CB**4.) 

VFB=VFB/E/AI 

RETURN 

ci  1 1 1 1  mi  1  mini 11 11  m  1 inn I V u  i  um/m  mi  um  i  m  mu  m  u 


a  1 1 1  u  1 1  u  1 11  / 11 1 1 1  m  11  / '  1 1  m  1 1  m  /  u  u  u  m  1  /  m  u  1 1 1 1  m  /  / 1  u  1  / 1 

FUNCTION  VSM(X,B#ALf E,AI , AM,DfG) 
C      DEFLECTION  FOR  REAM  WITH  A  SIMPLE  PENDULLM 

Q  *  if  *  *  *  *  if  if  it  *  «  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  * 

C****»  *  ************  -4  ****************************  ********************** 

BX=3*X 

AG=386. 

BL=B*AL 

C11=C0SH(BL)*C0S(BX)+C0S{BL)*C0SH(BX) 

C12=Cll/(C0SiBL)+C0SH(BL) ) 

C1=C12-1. 

GS=G*G 

Dll=i./GS-D/AG 

Dl=-AK/Dll 

02= E* A I *B*B*B/ (COS ( BL ) +C0SH ( BL ) ) 

D2=D2*386. 

C21=2.+2.*C0SH(BL)*C0S(BLJ 

C22=D2*C2i 

C23={SIN(RL)*SINH(BL) ) * ( COS ( BL )-C0SH( BL ) ) 

C2A=C23/(C0S(BL)+C0SH(8L) ) 

C25=SIN(BL)-SINH(BL)-C24 

C2=0l*C25+C22 

C3  1=  (  S  I N ( BL }  +  S I NH ( BL ) ) * ( COS  < BX } -COSH  (  BX ) ) 

C32=C31/(COS(BL)+COSH{BL) ) 

C 33=  S I N ! BX ) -S I NH ( BX ) -C 32 

C34=CCSHv3L)*SIN(BL)+COS(BL)*SINH(BL) 

C35=D2*C34 
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C36=i.-2.*C0SH(BL)*C0S{BL)/(C0S(BL)+C0SH(BL)) 

C37=C36*D1 
C3=(C37-i-C35)*C33 
VSK=(C3/C2+Cl)/(iJ**4.) 
VSM=VSM/E/AI 

RETURN 

C/  /  /  /  /  //// // // '// 7/ it  It I '// '// 7/ 7/ 7 7/ 7 '// 7/ 7 7/ 7 7 7 7 '// 7 7 7/ 7 7/ '// 7/ 7 // 7 7 7/ 7 7 V 


Zltttttt 

su 
s 


C  *  *  V  v  >> 

C      THI 

C   W  I  T 
c  *  *  *  *  * 


/// 

BRO 
*  x  -f 

SUB 


101 


C 
C 


H 

jl  -v. 
V  v 

d: 
di 
en 

FC 
AL 
CO 
DO 
DO 

AM 


CO 

DO 

DO 

SA 

CO 

K 

AN 

AN 

AN 

AN 

IF 

m= 

DO 
AN 
AN 
AN 


3 
25 


A  M 
*** 

KEN 

MEN 

MMO 

RMA 

SO 

HPU 

11  = 

1J= 

(Iv 

NT  I 


21  = 
2J  = 
M  {  I 
NT  I 
IS 
I  If 
i  If 
(K, 
IK, 
iK. 
K-l 
31  = 
(I. 
v  i  » 
(I. 


//// 
UTIN 

if  if  if  if 

ROUT 

ULTI 

„*«  j.  -v.  ^*. 
v  V  -.-  «P 

SI  ON 

SI  CM 

N  AM 

T  I  1 H 

NEED 

TES 

ifK 

l,K 

J)=0 

NUE 

1 » K 

IfK 

•  J)  = 

NUE 

THE 

1)=- 

2)  =  A 

K-l  J 

K)=- 

EQ.2 


////////// 

E  VM 
if  if  *  * 

I  ME 
PLE 


1PU.B, 

CALCUL. 
PENDUL 

:  <:  if  t=  if  if  if  if 


VF( 
AN( 
(  10) 
0,5E 
S  FU 
DEFL 


20),VP 
10, 10) 
,CL(10 
14. '3) 

NCTION 
ECTION 


//////////////////////// 
ALfEt'AIfGtOltVNP) 

if  if  if  if  if  if  if  *  *  *  *  *  *  *  <<  ♦  *  *  *  *  *  *  *  * 

ATES  THE  DEFLECTION  OF  A 
UM  ABSORBER. 

if  if  if  if  if  if  if  if  if  if  *  *  *  #  *  *  *  #  *  *  *  *  *  * 

P(  10) , VAB(20) 
tSAM(iO,10),R(10,10),N(l 

) , AR(IOO) ,K 

SUBPROGRAM  AAM(I) 

FOR  BEAM  WITH  MULTIPLE 


//////////////////// 


if  if  $  if  if  )(< 

CANTI 

0),L(1 


if  if  if  if  if 

LEVER 


if  if  if  it 

SEA 


if  if  if  if  v  V 
M 


i,:  if  if  *  if  if  if  if  if  if  it  if  if  *  <* 


0) 


PENDULUUM 


0.0 

ORDER  OF  T 
AAH(1)/CL( 
AM(2)/CL(2 
=AM{K)/CL( 
AM(K)/CL(K 
)G0  TO  25 


HE  PENDULUM 

i)-AAM(2)/CL(2) 

) 

K) 

) 


2,M 

I-1)=AAM{ I )/CL 
I)=-AAM(I )/CL( 
I+1)=AAM( I+l)/ 


(I) 

I)-AAMU  +  1)/CL(I  +  1) 

CL(I+1) 


CONTINU! 

CONTlNUi 

D04 

SAM 

CON 

DC  9 

D05 


CON 
CON 
D01 


I  =  1  ,  K 

(1,1 )=AMl I ) 

TINUE 

1  =  1,  K 

J=1,K 

, J)=-G*G*SAMU, J)-386.*AN( I, J) 

TINUE 

T  INUE 

01  =  1,  K 
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D06J=1?K 

6  AR{ (J-l>*K+I )=R( I, J) 
10  CONTINUE 

CALL  MINV(AR,K, D,L, N) 

H=0. 

D07I=1,K 

H=H+(AM(I) )*AR(I ) 

7  CONTINUE 
D1=-G*G*386.*AAK(1}*H/CL(1) 
BX=B*X 

AG=386. 

BL=B*AL 

CU=C0SH(BL)*C0S(BX)+C0S{BL)*C0SH{8X) 

C12=Cli/(COS(BL)+COSH(BL) ) 

Cl=Cl2-i. 

GS=G*G 

D2=E*AI*B*B*3/(C0S(BL)+C0SH(BL)  ) 

D2=D2*386. 

C2  1  =  2 .  +  2 . -COSH { BL  J  *COS ( BL ) 

C22=D2*C21 

C23= ( S IN( BL ) +SINHI BL ) )*( COS ( BL )-COSH ( BL ) ) 

C24=C23/(COS(BL)+COSH(BL) > 

J25=SIN(BLi-SINH(BL J-C24 

C2=0L*C25+C22 

C31={SIN(BL»+SINH(BL) ) * (COS { BX }-COSH( BX } ) 

C32=C31/{COS<BL)+COSH(BL)  ) 

C33*SIN(8X)-SINH<8X)-C32 

C34=C0SJ).iBL)*SI\,(3L)+C0S{BL)*SINH(BL) 

C35*02*C34 

C36=1.-2.*C0SH{3L)*C0S(BL)/(C0S(BL)+C0SH(BL) ) 

C37=C36*D1 
C3=C2/(C37+C35J 

IF .(X.LT.55.  5CO  TO  8 
WRITE(3»101)C8,D,C2,C37,C35 

8  CONTINUE 
C3=(C37*C35)*C33 

VNP=< (C3/C2)+Cl)/B**4./E/AI 

RETURN 
END 

:///// /////////////////////////////////////////////////////////////// 


c/// ////////////////////////////// /////////////////////////////////// 

FUNCTION  AAM{ I ) 

C#*«  $  ****** ********************#*************************************=* 

C   THIS  SUBROUTINE  CALCULATES 

C   A X  i I ) +AM ( I  +  L )  +  • . . +AK( K ) 

Q ***  a************************************************** *************** 

COMMON  AM (  1 0 ) , CL ( 1 0 ) , AR { 1 00  J , K 
AAM=0. 
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D01 1  1  =  1 ,K 
AAM=AAM+AK{ I  1) 
1  CONTINUE 
RETURN 
END 
C/  ///  /  / 7 I 7/  / 1 I  / I  ///  /  // 1 It It I II I I  It II I 7 '// 7 '// 7 7 7V 7 '// 7 7 7/ V 7 7 7/ 7/ 7 7/  / 7 7 7/ 

C/// /////////////////////// III  II llllt It  II 1 1  till  1 1 II II II II  tut  It  I  III  II 

c 

C   THIS  SUBROUTINE  WAS  TAKEN  FROM  IBM  SYSTEM/360  SCIENTIFIC  SUBROUTINE 

C   PACKAGE  (360A-CM-03X)  VERSION  II  PROGRAMED  S  MANUAL. 

C      

c 

C         SUBROUTINE  MINV 

C 

C  '      PURPOSE 

C  INVERT  A  MATRIX 

C 

C         USAGE 

C  CALL  MINV(A,N,D,L,M) 

C         DESCRIPTION  OF  PARAMETERS 

C  A  -  INPUT  MATRIX,  DESTROYED  IN  COMPUTATION  AND  REPLACED  BV 

C  RESULTANT  INVERSE. 

G  N  -  ORDER  OF  MATRIX  A 

G  0  -  RESULTANT  DETERMINANT 

G  L  -  WORK  VECTOR  OF  LENGTH  N 

C  M  -  WORK  VECTOR  OF  LENGTH  N 

C 

G         REMARKS 

G  MATRIX  A  MUST  BE  A  GENERAL  MATRIX 

C  . 

C         SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C  NONE 

C 

C  THE  STANDARD  GAUSS-JORDAN  METHOD  IS  USED.  THE  DETERMINANT 

C  IS  ALSO  CALCULATED.  A  DETERMINANT  OF  ZERO  INDICATES  THAi 

C  THE  MATRIX  IS  SINGULAR. 

C 

C  

c 


SUBROUTINE  M I NV ( A, N , D, L , M ) 

DIMENSI0NA(36) ,L(6) , M ( 6  J 
C 
C        

c 

C         IF  A  DOUBLE  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED,  THE 
C         G  IN  COLUMN  1  SHOULD  BE  REMOVED  FROM  THE  D0U3LE  PRECISION 
C         STATEMENT  V.HICH  FOLLOWS. 
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C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 
c 
c 


D0U3LE  PRECISION  A, D, BIGA v HOLD 

THE  C  MUST  ALSO  BE  REMOVED  FROM  DOUBLE  PRECISION  STATEMENTS 
APPEARING  IN  OTHER  ROUTINES  USED  IN  CONJUNCTION  WITH  THIS 
ROUTINE. 

THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
CONTAIN  DOUBLE  PRECISION  FORTRAN  FUNCTIONS.   ABS  IN  STATEMENT 
10  MUST  BE  CHANGED  TO  DABS. 


SEARCH  FOR  LARGEST  ELEMENT 

D=1.0 

NK=-N 

DO    SO    K=1,N 

NK  =  NK-rN 

L(K)=K 

M  {  K  )  =  K 

KK=NK+K 

biga=a;kk) 

DO  20  J=K,N 
IZ=N*(J-1) 
DO  20  I=K.,N 
IJ=IZ+I 

io  ;f(  abs(biga)-  abs(aiij)))  15,20,20 

15  3IGA=A<IJ) 
L(K)=I 

K  (  K  J  =  J 

20  CONTINUE 

INTERCHANGE  ROWS 

J=L(K) 

IF(J-X)  35,35,25 

25  KI=K-N 

DO  30  1=1, N 

k.i=ki+n 

HOLD=-A(KI ) 
JI=KI-K+J 

A(KI)=A(JI) 
30  A(JI )  =  HOLD 


C 

C 


INTERCHANGE  COLUMNS 


35  i=m;kj 

Ir(I-K)  45,45,38 

3  0     r\ if  j.  f  y    ■*  ' 
0   J  r>  =  r.  V  k    1  —  i  / 

DO  4  0  J=1,N 

jk=n;<+j 


■> 


6o 


JI=JP+J 
HCLD=-A(JK) 
A(JK)=A(JI) 
40  A(JI  )  =H0L0 

C  -    T, 

C        DIVIDE  COLUMN  BY  MINUS  PIVOT  (VALUE  OF  PIVOT  ELEKtM  *S 

C         CONTAINED  IN  BIGA) 

c 

45  IF(BIGA)  48,46,48 

46  D  =  0.0 
RETURN 

43  DO  5  5  1  =  1  ,N 

IF(I-K)  50,55,50 
5C  IK=NK+I 

At [K)=Al IK)/(-8IGA) 
55  CONTINUE 

C 

C         REDUCE  MATRIX 

C 

DO  65  1  =  1, N 

IK=NK+I 

HOLD=A( IK) 

IJ=I-N 

DO  65  J=1,N 

IJ=IJvN 

IFII-X)  60,65,60 
60  IFIJ-K)  62,65,62 
62  KJ=IJ-I+K 

A( IJ)=HOLD*A{KJ)+A( IJ) 

AUJ)=AUK>*A(KJ)  +  A(IJ) 
6  5  CONTINUE 
C 
C         DIVIDE  ROW  BY  PIVOT 

C 

i<J  =  K-N 
DO    75    J=1,N 
KJ=KJ+N 

IF(J-K)  70,75,70 
70    AC<J;=AC<J)/BIGA 

75  CONTINUE  . 

C 
C         PRODUCT  OF  PIVOTS 

C 

D=D*8IGA 

C 

C         REPLACE  PIVOT  BY  RECIPROCAL 

C 

A(KK)=1.0/bIGA 

80  CONTINUE 

C 

C         FINAL  ROW  AND  COLUMN  INTERCHANGE 


61 


K  =  N 

100  K=  (;<-!> 

IRK)     150,150,105 
105    I=UK) 

IFU-K)     120,120,108 
108    JQ=N*(K-1) 

JR=N*( 1-1) 

DO    110    J  =  l,fl 

JK=JQ+J 

HOLO=A(JK) 

JI=JR+J 

A( JK)=-A( J  I  ) 
110    AUI)    =H0LD 
120    J=K(K) 

IF  ( J-IO     100,100,125 
125    KI=K-N 

DO    130    1=1, N 

KI=KI+N 

HOLD=A(KI ) 

JI=KI-K+J 

A(KI }=-A(JI ) 
130    AUI  )    =HOLO 

GO    TO    100 
150    RETURN 

^.i  1 1 1  nfi  1 1 1  n  1 1  n  j  i  u  i  n  i  n  1 1 1 1 1  n  n  n  n  1 1 1 1  n  n  m  1 1  / 1 1  n  1 1  n  1 1 1 1 1 1 


PROCESSING    ENDED    AT    EOD 
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